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Abstract 

X-ray phase-contrast imaging sliows improved soft-tissue contrast compared to standard absorption-based X-ray imaging. 
Especially the grating-based method seems to be one promising candidate for clinical implementation due to its 
extendibility to standard laboratory X-ray sources. Therefore the purpose of our study was to evaluate the potential of 
grating-based phase-contrast computed tomography in combination with a novel bi-lateral denoising method for imaging 
of focal liver lesions in an ex vivo feasibility study. Our study shows that grating-based phase-contrast CT (PCCT) significantly 
increases the soft-tissue contrast in the ex vivo liver specimens. Combining the information of both signals - absorption and 
phase-contrast - the bi-lateral filtering leads to an improvement of lesion detectability and higher contrast-to-noise ratios. 
The normal and the pathological tissue can be clearly delineated and even internal structures of the pathological tissue can 
be visualized, being invisible in the absorption-based CT alone. Histopathology confirmed the presence of the 
corresponding findings in the analyzed tissue. The results give strong evidence for a sufficiently high contrast for different 
liver lesions using non-contrast-enhanced PCCT. Thus, ex vivo imaging of liver lesions is possible with a polychromatic X-ray 
source and at a spatial resolution of ~100 ^m. The post-processing with the novel bi-lateral denoising method improves 
the image quality by combining the information from the absorption and the phase-contrast images. 
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Introduction 

Focal liver lesions are a common finding in ultrasound and 
computed tomography (CT) examinations often presenting 
without clinical symptoms [1-3]. Whether a hepatic lesion 
requires treatment, follow-up imaging, or does not need further 
attention depends on the dignity (benign vs. malignant) and entity 
of the finding [4] . Furthermore, in oncology the liver is the most 
common site of metastasis from gastrointestinal tumors. Accurate 
assessment of tumor spread and treatment response are important 
for therapy planning [5] . 

Non-enhanced CT is very limited for the detection and 
characterization of focal hepatic lesions due to its low soft-tissue 
contrast. The use of iodine contrast agents greatly improves 
detectability of lesions; however, the characterization especially of 
subcentimetric findings remains challenging [6]. Also, in some 
patients iodine-based contrast media can cause anaphylactic 
reactions, acute renal failure, or hyperthyroidism [7-9]. Drug 
interactions with metformin, an oral antidiabetic agent, have been 
described [10]. The gold standard for the detection and 



characterization of focal hepatic lesions remains magnetic 
resonance imaging (MRI) with the application of diffusion 
weighted imaging (DWI) and liver-specific contrast agents 
[11,12]. Contraindications for MRI include cardiac pacemakers 
and patients suffering from claustrophobia. Additionally, long 
scanning times can be a limitation due to motion artifacts or 
limited patient tolerance, as well as the limited spatial resolution 
compared to CT. 

Conventional multi-slice CT uses the attenuation of x-rays by 
photo-electric absorption and Compton scattering processes in the 
examined object as the sole source of contrast. Therefore, 
structures containing elements with high Z number, like calcium 
in bones, give high contrast whereas soft tissues such as muscle, 
parenchymatous organs and adipose tissue show rather low 
contrast. On the other hand, phase-contrast computed tomogra- 
phy (PCCT) exploits the physical effect of the phase shift x-rays 
experience when passing through an object. 

Different phase-contrast methods have been developed over the 
last decades [13]. All these methods can be divided into three main 
groups: (1.) the direct methods measuring the phase shift (Bonse- 
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Hart interferometer) [14,15], (2.) the propagation based methods 
measuring the 2nd derivative of the phase shift [16], and (3.) the 
differential methods measuring the 1 st derivative of the phase shift 
by using e.g. the analyzer-based [17,18] or the grating-based 
method [19,20]. Only the grating-based phase-contrast computed 
tomography (PCCT) has been successfully expanded to standard, 
polychromatic X-ray tubes [21]. This method provides a three- 
dimensional map of the attenuation coefficients and the electron 
density distribution of the specimen at the same time. Several 
recent articles have shown the potential of the grating-based 
PCCT for clinical applications demonstrating a remarkable 
improvement of soft-tissue contrast without the use of contrast 
media [22-25]. 

The purpose of our experimental proof-of-principle study was to 
evaluate the potential of grating-based PCCT for the detection 
and characterization of focal lesions in ex-vivo human liver 
specimens using a standard laboratory x-ray tube setup. Further- 
more, we aimed to apply a novel post-processing denoising 
method on the tomographic data to demonstrate the possible 
quality improvement of the images. 

Materials and Methods 

Grating interferometry 

The principle of grating-based PCCT and its projection 
acquisition is explained in detail in [19-21]. The x-rays are 
attenuated and refracted when they pass through the object. The 
refraction causes a change in the direction of the x-ray path, which 
only can be measured indrrecdy using for example an interfer- 
ometer as shown in Figure 1. This method provides the phase- 
contrast and the absorption-contrast images simultaneously. The 
device used with coherent synchrotron sources consists of two x- 
ray optical gratings [26,27]: the phase grating produces an 
interference pattern in discrete distances downstream the x-ray 
beam, and the analyzer grating is used to detect the pattern of 
several micrometer periods with a standard x-ray imaging detector 
by scanning the analyzer grating with respect to the phase grating 
[21]. The method can be transferred to standard laboratory 
sources by adding an additional source grating, which increases 
the coherence of the X-ray source. 

Experimental setup 

For the experimental setup a three grating interferometer has 
been used in combination with an ENRAF Nonius rotating anode 
X-ray tube and a single photon counting detector (Pilatus II, 
Dectris, Baden, Switzerland) with a pixel size of 172x172 (im 
and a silicon sensor thickness of 450 |J,m. Gratings with periods of 
5.4 |Xm (Microworks, Karlsruhe, Germany) in a symmetric setup 
have been chosen: Both inter-grating distances have been set to 
875 mm. The phase grating made of Nickel had a height of 8 |J.m 
and introduced a phase-shift of IT at the design energy of 
22.8 keV. The golden source and analyzer grating structures were 
50 |J,m high. The 6'* Talbot order for the design energy of 23 keV 
was chosen and the mean visibility of the interferometer reached 
13%. The X-ray filtration consisted of the following components: 

• Beryllium window of the X-ray tube with a thickness of 
0.25 mm 

• Silicon wafers of source, phase and analyzer grating with a 
thickness 0.5 mm each 

• Silicon disk with a thickness of 1 mm 

• Approximately 2.25 m of air. 



Beryllium window and silicon wafers of the gratings are not 
chosen based on advanced simulations but are parts of the 
individual components of the set-up. The silicon disk has been 
introduced to further reduce low energy photons without relevant 
reduction of photon flux. 

The target material of the X-ray tube was Molybdenum and a 
tube voltage of 35 kVp has been applied in all tomography scans. 
The tube current has been set to 70 niA resulting in a total tube 
power of 2.45 kW. Altogether, without an object in the beam, an 
average count of 630 photons per second has been collected in 
each pixel of the detector. The actual size of the focal spot is 
3x0.3 mm^, the effective source size is only 0.3x0.3 mm^ due to 
an inclination of the anode surface with respect to the optical axis 
of the beam. 

The whole interferometer has been installed 60 cm away from 
the source, each sample was placed 70 mm in front of the phase 
grating and the detector 40 mm behind the analyzer grating, 
which gives a sample magnification of 1.72 and an effective pixel 
size of 0.1 xO.l mm^. 

Liver samples 

AH samples were acquired through the Institute of Padiology, 
Technische Universitat Mtinchen. They were initially resected for 
clinical reasons independently of this study. The patients provided 
their written consent to the use of tissue samples for educational 
and scientific purposes in general. The use of these specimens and 
the way how written consent was obtained has been fully approved 
by the ethics committee at the Faculty of Medicine of the 
Technische Universitat Mtinchen. 

After excision a selected representative tissue sample was excised 
and put in a 50 ml plastic container in a 4% neutral-buffered 
formaldehyde solution. The following five different liver specimens 
with lesions have been chosen: 1.) a metastasis of a low-grade 
adenocarcinoma in a steatotic liver, 2.) a cholangiocellular 
carcinoma in a liver with macrosteatosis and pilosis, 3.) a 
hepatocellular carcinoma in a cirrhotic liver, 4.) a subcapsular 
hematoma and 5.) a metastasis of a mucinous adenocarcinoma of 
the colon. Histological sections of the specimens were prepared 
with a standard hematoxylin and eosion (H&E) staining. The slice 
thickness was 5 micrometers. Histology and PCCT slices were 
manually correlated using prominent features for the correct 
orientation. 

Image acquisition 

The liver samples were put into cylindrical plastic containers of 
30 mm diameter that were scanned in air. 1200 projections over 
360° were recorded per tomography scan. Every 20 projections 5 
reference projections without the sample were recorded for 
background correction. Each projection was obtained by the 
phase stepping method [19] using eleven images with 5 seconds 
exposure time acquired over one period of the source grating. The 
total scan time was approximately 25 h. 

Image reconstruction and post processing 

The beam divergence is approximately 1.1° and allows to 
perform a standard parallel beam filtered back-projection using a 
Ram-Lak filter to reconstruct the absorption-contrast data and a 
Hubert filter to reconstruct the phase-contrast data [21]. The 
reconstructed slice thickness of 100 |xm is given by the effective 
pixel size of the detector, and the vertical number of detector 
pixels limits the number of slices to 195. Joined bilateral filtering as 
described in [28] was used as post processing in order to reduce 
the noise level. This method exploits the fact that the absorption 
and phase-contrast images are perfectly registered with respect to 
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Figure 1. Schematic and a photograpKi of the grating-based imaging setup. A) Schematic top view of the grating-based imaging setup. The 
X-ray tube is a rotating anode with a Molybdenum target working with max. 60 kVp and 80 mA. The grating-interferometer consists of a source 
grating (to increase the coherence), a phase grating, an analyzer grating and a detector. The specimen is rotated in this setup. B) Photograph of the 
grating-based PCCT setup using conventional X-ray sources consisting of three gratings arranged in a symmetrical setup: (a) source grating marked 
yellow, (b) phase grating marked purple and (c) analyzer grating marked red, (d) the specimen mounted to a hanging rotation stage and (e) the 
detector. A red arrow marks the X-ray beam path. 
doi:10.1371/journal.pone.0083369.g001 




Figure 2. Comparison of unfiltered and filtered tomographic slices of a metastasis of a low-grade adenocarcinoma in a steatotic 
liver. Axial slices of the conventional absorption-based (A, D) and the phase-contrast (B, E) tomography windowed using the same window level and 
window width. The histograms of both unfiltered (C) and filtered (F) tomographic datasets demonstrate the filtering result and the red dashed 
markers show the window level. The areas I (red, surrounding formalin) and II (blue, high-contrast tumor in (B) mark the regions of 30x30 pixels, 
which were averaged for CNR calculation. 
doi:1 0.1 371/journal.pone.0083369.g002 
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Table 1. Comparison of CNRs for unfiltered and filtered tomographic slices of a metastasis of a low-grade adenocarcinoma in a 
steatotic liver shown in Figure 2. 





absorption 


absorption phase-contrast phase-contrast 


unfiltered 


filtered unfiltered filtered 


CNR between ROI 1 and II 0.39 


1.24 5.34 8.63 


For the calculation of the CNRs the areas 1 (red, surrounding formalin) and 


(blue, high-contrast tumor) marked in Figure 2B were used, which correspond to 30x30 



pixels. 

doi:l 0.1 371 /journal.pone.0083369.t001 



each other by performing edge-preserving filtering where edges 
are detected simultaneously in both images. This method 
overcomes at least partially the problem of standard edge- 
preserving filtering where edges with a contrast-to-noise-ratio 
(CNR) less than 1 are typically not detected as edges and are thus 
smoothed. Due to the simultaneous detection of edges in both 
images, edge detection becomes more reliable and thus, a mutual 
benefit is created. For the liver samples here, the noise levels, 
which are free parameters in the joined bilateral filter, were 
estimated as the local standard deviation of reconstructed image 
values in homogeneous regions. The filter was applied using a 1 3' 
neighborhood. The contrast-to-noise ratios (CNR) of two different 
regions (30x30 pixels each) of one sample relative to the 
surrounding formalin were calculated in absorption and phase- 
contrast tomogram before and after the denoising. The regions 
used for the CNR calculation are marked in Figure 2B The 
following equation was used to determine the CNRs: 



CNR- 



'Sh\ 



OS 



(1) 



where Sa and Sj represent the averaged measured signals in the 
region a and A, respectively, and as is the error of the averaged 
signal. Here we set as = a,,, which is the standard deviation of the 
averaged signal Sa outside the specimen, to avoid higher noise 
values when averaging over inhomogeneous tissue. Please note 
that this implies the assumption that the noise is equally distributed 
over the whole image and the standard deviation cr„ of the 
homogenous region outside the specimen represents the noise level 
in the image. 

Results 

Figure 2 compares the unfiltered and the filtered tomographic 
slices of a metastasis of a low-grade adenocarcinoma in a steatotic 



absorption 



phase contrast 



c 
O 




liver tissue 



Figure 3. Filtered tomographic slices of a metastasis of a low-grade adenocarcinoma in a steatotic liver filtered with the 3D 
bilateral filter. Conventional absorption-based (A, D) and phase-contrast (B, E) CT images of the liver specimen in two different planes. H&E stain (C) 
corresponds to image A and B, respectively. (F) shows the histogram of the filtered absorption image (top) and the filtered phase-contrast image 
(bottom) with red dashed lines marking the window level in the shown slices. 
doi:10.1371/journal.pone.0083369.g003 
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Figure 4. Tomographic slices of a cKioiangiocelluiar carcinoma in a liver witKi macrosteatosis and pilosis filtered with the 3D 
bilateral filter. Conventional absorption-based (A, D) and phase-contrast (B, E) CT images of the liver specimen in two different planes. H&E stain (C) 
corresponds to image A and B, respectively. (F) shows the histogram of the filtered absorption image (top) and the filtered phase-contrast image 
(bottom) with red dashed lines marking the window level in the shown slices. The absorption images (A, D) profit strongly from the filtering, but the 
information is still complementary. 
doi:1 0.1 371 /journal.pone.0083369.g004 



liver. For this comparison the same window level and window 
width was used for the slices of both contrasts, which is indicated 
by red dashed lines in the histograms of the unfiltered (C) and of 
the filtered shces (F). The unfiltered absorption images (A) show a 
very low soft-tissue contrast and a high noise level. The areas with 
steatotic liver tissue show low density compared to the surrounding 
formalin, while the infiltrating tumor tissue in the center has a 
higher density. In the phase-contrast images (B, E) the steatotic 
liver tissue has also a low signal, but the tumor tissue shows 
intermediate to high signal with more details compared to the 
absorption images. Both signals profit from the post processing 
with the bilateral filter, but especially the absorption-contrast 
image contains many structural details inside the tumor after the 
filtering. For the remaining samples only the filtered images will be 
shown. To quantify the image-quality improvement the CNR of 
one region inside the tumour compared to the formalin was 
calculated for both signals. Note, that the two tomograms are 
intrinsically registered, thus, the marked regions-of-interest (ROIs) 
I (red, surrounding formalin) and II (blue, high-contrast tumor 
tissue) in Figure 2B correspond to 30 x30 pixels and are exactly the 
same for the slices shown in Figures 2A, 2D, and 2E. These ROIs 
were averaged to determine the CNRs as described in the methods 
and materials section, where the signal S,, was chosen from the 
ROI I and the signal Sb from the ROI II. The resulting CNRs are 
given in Table 1. 

In Figure 3 only the filtered tomographic slices of the same 
metastasis as shown in Figxire 2 are depicted. The conventional 
absorption based (3A, 3D) and phase-contrast (3B, 3E) CT images 
of the liver specimen are shown in two different planes axial and 
frontal. The H&E stained histopathology image (3C) shows a 



corresponding region of the specimen like in 3A and 3B, 
respectively. The histograms (3F) of both filtered signals are 
presented to indicate the chflFerent windowing of the images, which 
is set to a visual optimum for each signal. The absorption images 
(3A, 3D) profit mostiy from the post-processing in the higher- 
density area of the infiltrating tumor tissue in the center, while the 
areas with steatotic liver tissue with low density lead to an CNR 
improvement of this region in the phase-contrast image. 

Figure 4 shows filtered tomographic slices of a cholangiocellular 
carcinoma in a liver with macrosteatosis and pilosis. (4A, 4D) 
present the conventional absorption-based and (4B, 4E) the phase- 
contrast CT slices of the liver specimen in two different planes. In 
(4C) the H&E stained slice corresponding to the regions of images 
4A and 4B is shown. Even after the bilateral filtering the soft-tissue 
contrast in absorption images (4A and 4D) is low in comparison 
with the phase-contrast images (4B, 4E), which show low signal for 
the tumor tissue in the center of the slice with high signal in 
necrotic/hemorrhagic areas compared to the normal liver tissue in 
the upper left and lower right parts (4B) with intermediate signal. 
Here, the windowing is chosen to a visual optimum of each signal 
and is indicated in the histograms (4F). 

The filtered tomographic slices of a hepatocellular carcinoma in 
a cirrhotic liver are shown in Figure 5. The patient had received a 
transcatheter arterial chemoembolisation (TACE) prior to exci- 
sion. In the absorption images (5A, 5D) the HCC nodules (marked 
with red arrows in 5A) in the right center are visible due to the 
high density of the retained Lipiodol. Lipiodol is a composition of 
iodine with poppyseed oil, nowadays mainly used for chemoem- 
bolizatioii procedures where a chemotherapeutic agent is mixed 
with the oily liquid. Chemoembolization interrupts the tumors 
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Figure 5. Filtered tomographic slices of a hepatocellular carcinoma in a cirrhotic liver. The patient had received transcatheter arterial 
chemoembolisation (TACE). Conventional absorption-based (A, D) and phase-contrast (B, E) CT images of the liver specimen in two different planes. 
H&E stain (C) corresponds to image A and B, respectively, and (F) shows the histogram of the absorption (top) and the phase-contrast image (bottom) 
with red dashed lines marking the window level. In the absorption images (A, D) the HCC nodules in the right center (marked with red arrows in A) 
are visible due to the high density of the retained Lipiodol. In the phase-contrast images (B, D) the fibrous septa are visible as high signal bands 
whereas the liver tissue and the HCC nodules show intermediate signal. Some areas with increased Lipiodol retention show low signal. 
doi:1 0.1 371 /journai.pone.0083369.g005 



blood supply and allows high local concentrations of the drug for a 
longer period of time. In the phase-contrast images (5B, 5D) the 
fibrous septa are visible as high signal bands whereas the liver 
tissue and the HCC nodules show intermediate signal. Some areas 
with increased Lipiodol retention show low signal. Again the 
windowing is set to the visual optimum of each signal, which is 
marked in the histograms (5F). 

Figure 6 presents the tomographic slices of a subcapsular 
hematoma after bilateral filtering. The contrast between the liver 
tissue in the upper right part of the images and the subcapsular 
hematoma in the lower left part is considerably lower in the 
absorption images (6A and 6D) compared to the phase-contrast 
images (6B and 6E), where the hematoma shows a high signal. 

The last specimen was a metastasis of a mucinous adenocar- 
cinoma of the colon. The corresponding filtered tomographic 
shces are shown in Figure 7. The phase-contrast images (7B, 7E) 
show a significantly higher soft-tissue contrast with low signal in 
the mucinous areas, intermediate signal in the liver tissue and 
higher signal for the tumor tissue and necrotic/hemorrhagic areas. 
In the absorption images (7 A, 7C) only the mucinous parts of the 
metastasis are clearly visible as areas with low density. The other 
structures in the absorption images are barely visible since the 
contrast after filtering is in the same order of magnitude as the 
remaining low-frequency artifacts. Nevertheless, sharp edges are 
visible, which have been rehably detected in the phase-contrast 
image. 



Discussion 

Our results clearly demonstrate that the tumor boundary and 
the soft-tissue in general can be visualized with a lab-based 
imaging setup. The effective pixel sizes of 100 |J,m are one 
magnitude larger than the previously reported experimental results 
on human liver tissue [29] and on other tissue types gained at 
brilliant synchrotron sources [25,30] . This is not yet in the range of 
clinically relevant spatial resolution (i.e. 500 |J.m instead of 100 |im 
in our study) required to meet the stringent dose requirements for 
in-vivo applications [31], but it is the first experimental 
demonstration of the stiU present information gain when using 
PCCT in combination with polychromatic radiation and low 
spatial resolution for ex vivo human liver samples. 

By the use of joined bilateral filtering, the CNR in both images 
was improved substantially. In the present study, most of the 
originally reconstructed phase-contrast images showed a much 
better CNR than the absorption-contrast images. Therefore, there 
was less mutual benefit but more a unidirectional one, viz. the 
absorption-contrast image gained much more than the phase- 
contrast image: Edges, which were completely invisible in the 
original FBP images were made visible thanks to their visibility in 
the phase-contrast image. However, we also found some high- 
contrast regions in the absorption-contrast image having a better 
CNR and in this region, the difiFerential phase-contrast image 
profited from the joined processing. With regard to the different 
visible structures inside the liver specimens, especially the 
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Figure 6. Filtered tomograpKiic slices of a subcapsular hematoma. Conventional absorption- based (A, D) and phase-contrast (6, E) CT images 
of the liver specimen in two different planes. H&E stain (C) corresponds to image A and B, respectively, and the histogram (F) of the absorption (top) 
and the phase-contrast image (bottom) show the window level as red dashed lines. In the absorption images (A, D) the contrast between the liver 
tissue in the upper right part of the image (A, D) and the subcapsular hematoma in the lower left part is considerably lower compared to the phase- 
contrast image (B, E) where the hematoma shows a high signal. 
doi:1 0.1 371 /journal.pone.0083369.g006 



combination of botli information has the potential to improve tlie 
diagnostics without any use of contrast agents in future. 

The low-frequency noise (rings), which can be observed in the 
attenuation images, is a combination of a beam-hardening effect 
and remaining fringes in the background, which are not removed 
by the normalization with a flat-field image. We did not apply any 
corrections for these effects to our data. 

Our system has some limitations, which need to be discussed in 
more detail. One of the limitations of our study concerns the 
formalin fixation of the specimens. Especially this fact might 
change the contrast for different tissue types compared to in-vivo 
imaging of the tissue surrounded by circulating blood. Further 
studies are required to investigate the influence of formalin fixation 
on the contrast compared to fresh tissue in both, absorption and 
phase-contrast data. 

The second important limitation is related to the fact that our 
laboratory setup is still in an experimental state. This means that it 
was not optimized for fast measurements, which makes it 
comparable to MRI with a microscopy coil in terms of scanning 
times, but allowing better spatial resolution. Optimization of other 
parameters, which are important for the phase-contrast imaging 
(as for example filtering of the spectrum or setup geometry) are 
subject of ongoing studies, and are non-trivial in combination with 
a polychromatic source of a limited flux. These parameters 
strongly influence the quality of the interferometer and its 
sensitivity, but at the same time can reduce the number of 
photons, which downgrades the image quality. Thus, our 
measurements do also not meet any dose requirements for in- 
vivo applications. 



The apphed dose can be roughly estimated from the number of 
the photons recorded by the detector. Approximately 600 photons 
were collected per image, per second and detector pixel. Taking 
into account the absorption of photons by the analyzer grating 
(50%), the quantum efficiency of the PHatus detector (8% @ 
35 keV and 25% @ 22 keV) and the effective source spectrum, we 
can estimate that the total apphed dose to the sample was in the 
order of several Gy for a fuU tomography scan. Due to a non- 
optimized setup only a small fraction of the total applied dose 
really contributed to the signal formation and detection. Never- 
theless, the comparison of the presented images of both signals - 
absorption and phase-contrast - is fair, since both were obtained 
simultaneously and at exact the same dose. 

Various optimization steps apphed to our system can bring the 
applied dose per scan into the range of conventional micro CT 
systems (e.g. smaU-animal micro CT systems), which have 
comparably high spatial resolution and apply several hundreds 
of mOy for a full tomography scan. 

Optimization can be done from the detector side improving the 
detection efficiency of X-rays in the used energy range. The 
relatively low detector efficiency of our experimental system results 
in low photon statistic, which especially becomes apparent in the 
absorption images, where the image quality is worse than expected 
from measurements with an applied dose of several Gy. Modern 
detectors used in clinical systems can reach a quantum efficiency of 
up to 60% [32], which brings us already into the ball park of 
several hundreds mGy used in small animal scanning systems. A 
doubling of the visibility allows a dose reduction to a quarter at 
constant image quality [33]. Thus, raising visibility from currendy 
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Figure 7. Filtered tomograpKiic slices of a liver metastasis of a mucinous adenocarcinoma of the colon. Conventional absorption-based 
(A, D) and phase-contrast (B, E) CT images of the liver specimen in two different planes. H&E stain (C) corresponds to image A and B, respectively, and 
(F) shows the histograms of both signals with red dashed lines marking the window level. In the absorption images (A, C) only the mucinous parts of 
the metastasis are visible as areas with low density. The phase-contrast images (B, E) show a significantly higher soft-tissue contrast with low signal in 
the mucinous areas, intermediate signal in the liver tissue and higher signal for the tumor tissue and necrotic/hemorrhagic areas. 
doi:1 0.1 371 /journal.pone.0083369.g007 



10% to 30% by optimization of gratings and setup geometry 
would reduce dose to 1/9 of the current value. 

Advanced acquisition protocols can also reduce the acquisition 
time and the applied dose, which has been successfully demon- 
strated for grating-based imaging either theoretically or experi- 
mentally at synchrotron sources (low-dose grating-based phase- 
contrast imaging [34] and advanced stepping methods [35]). We 
also anticipate further possibilities to maintain the current CNRs 
at reduced dose and thus at reduced scan time by using iterative 
reconstruction techniques. In our experimental study the acqui- 
sition parameters like number of projections, number of phase 
steps, and acquisition time were arbitrarily chosen. Optimizing 
these parameters can significantly reduce the dose, e.g. utilizing a 
single-shot technique wiU allow a reduction of dose and acquisition 
time by a factor of 1 1 (as we have used 1 1 steps in the phase- 
stepping scan in our study) [36] . Furthermore, recently published 
results demonstrate that novel reconstruction methods can reduce 
dose and scanning time in phase-contrast CT of breast tissue by up 
to 74% [37]. 

The field-of-view of the grating-based imaging method is 
currently limited by the size of the available grating structures. 
Different efforts were made to overcome this limitation, by e.g. 
bending the wafer of large-area grating structures to account for 



the cone beam used for large field-of-views [38]. Recently, we 
could demonstrate that the grating-based method is able to work 
at high X-ray energies (up to 82 keV [39]), which brings it already 
in the range of clinically relevant X-ray energies. Further 
investigations are required to demonstrate the ability of the 
imaging method to image large objects in the range of the human 
body (20-40 cm diameter), which might lead to a degradation of 
image quality due to the non-negligible amount of scattering by 
large objects. 

In conclusion, with our work we demonstrate the feasibility to 
assess the potential of phase-contrast imaging for different kind of 
liver diseases on ex vivo samples using laboratory-based X-ray 
generators. Our method provides a three-dimensional map of tiie 
attenuation coefficients and the electron density distribution of the 
specimen at the same time and the joined filtering combines the 
information to improve the image quality. 
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